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1. INTRODUCTION 

Rubeosis iridis is an eyes disease induced from contemporary formation of abnormal blood vessels 
via surface of the iris through neovascularization [1] as shown in Figure 1. Normally, this disease is followed 
by neovascular glaucoma as it is one of the most common postoperative complications following vitrectomy 
for diabetic patient [2]. Typically, the early phase of this disease can be detected through ophthalmological 
examinations with aids of technologically-advanced instruments such as the fundus photography [3] and 
fluorescein angiography [4]. The screening test involve manual assessment by ophthalmology using high 
technology instrument is complicated, costly and time consuming. Therefore, the pre-screening rubeosis 
iridis system is highly needed as a tool prior a proper ophthalmologic consultation. 

This manual assessment by ophthalmology to examine the rubeosis iridis disease is really 
complicated and need more time to get the result. There is a simple, modern and non-intrusive pre-screening 
method to detect this disease through iris image processing system. Generally, this approach developed using 
several stages include image enhancement process, segmentation, feature extraction and classification [5]. 
However, this approach required high quality image in order to get the high accuracy results. In addition, 
there are various type of camera that can captured the retina image resulting multiresolution iris image. 
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Nowadays, current smartphone is not only used as a tool to communicate with others but also built in with camera. 
Therefore, user can easily use this smart phone camera to capture the iris image and doing the pre-screening of this 
disease before getting further ophthalmology consultation. However, the common problems that occured in the 
iris image processing system is the quality of image [6]. The image contrast can be low or high depending on 
environment factor such as environment, type of camera, surrounding condition and application used [7]. Also, 
if the image has poor quality, the amount of information required from the iris image may be reduced. 

Besides, image enhancement process is one of the crucial considerations in the iris image processing 
system. For last decades, many image enhancement methods have been studied and one of the popular 
techniques is enhancement process based on contrast. Basically, the philosophy of the contrast enhancement 
process i.e., to enhance the quality of image by enlarging the gray level image. The common of contrast 
enhancement process used is based on histogram equalization (HE) method [8]-[11] due to the simplicity and 
highly effectiveness. However, the major drawback of this method is losing the originality of image, loss of image 
information, over enhancement of brightness as well as the contrast and amplified the noise from the original 
image. Many attempts done by researchers to reduce the HE drawback by introducing several methods such as 
adaptive histogram equalization (AHE) [12], contrast limited adaptive histogram equalization (CLAHE) [13], 
brightness preserving bi-histogram equalization [14], sub-image histogram equalization method [15], recursive 
mean separate histogram equalization [16] and bi-histogram equalization [17]. 

Prior research in rubeosis iridis detection concentrated on image enhancement as a whole, without 
taking into account the images’ multi-contrast [18]-[22]. The process leads to over-enhancement or downgrade 
the quality of the image. Thus, this study proposed pre-enhancement process with the two well know contrast 
enhancement methods which is HE and CLAHE. Basically, the HE method will enhance the entire contrast 
image by enhancing the whole image histogram into the uniform histogram [8]. Meanwhile the CLAHE method 
modifies the contrast by each of small windows value based on some characteristics and it is applied to improve 
the contrast and minimize the noise that presence in the image [23]. Therefore, the main objective of this paper 
is to develop an algorithm to modify the contrast image in enhancement process by introducing the contrast 
modification technique for the pre-enhancement process, in order to detect rubeosis iridis disease. 


Figure 1. Iris with rubeosis iridis disease 


2. RESEARCH METHOD 
This study was designed using a numerical technique and quantitative analysis. Generally, this study 
consists of four stages as illustrated in Figure 2. The details of each stage are outlined. 


Dataset: Download 
rubeosis iridis images >] 
from online sources 


Image contrast Performance 
g — Enhancement — 
measurement measurement 


Figure 2. Block diagram for stages involved in the enhancement process 


2.1. Dataset 

The dataset was downloaded from publicly Google search engine [24]-[29]. Appropriated inputs 
then comprised 38 images of iris with rubeosis iridis disease from various web-based sources with 
multiresolution images range from (183x275) up to (1358x2048) pixels. For this study, images with an 
obvious abnormal blood vessel appeared were chosen as a dataset. Thus, it will lead to a high-performance 
measurement. The dataset was classified into three groups; low, medium, and high resolution. The classification 
is based on spatial measurements pixels per inch (ppi), width (w), and height (h). The (1) defines the group 
for each of the conditions of the resulted image, I: 
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300 ppi < w < 500 ppi 
low; I(w,h) = | and 
h « 1000 ppi 


w « 500 ppi 
medium; I(w, h) = | and 
500 ppi < h < 500 ppi 


w > 600 ppi 
high; I(w,h) = | and (1) 
h = 1000 ppi 


2.2. Image contrast measurement 

Figure 3 depicts the image contrast measuring procedure. The effectiveness of the luminance and 
arithmetic mean model was investigated to identify the appropriate model for the dataset. Luminance refers 
to the brightness of light emitted or reflected from a surface. The (2) shows a mathematical model for 
luminance, Y; : 


Y, = 0.299R + 0.587G + 0.114B (2) 


Where R, G, and B denote pixel values at the the respective red, green, and blue channels. Notably, 
parameters R and G outweigh parameter B [1] due to the human eye is more sensitive to the green and red 
colours. Meanwhile, arithmetic mean model which depicts the average luminance, Y of all points on a surface 
is outlined as in (3): 


_ (R+G+B) 
= 3 


Y (3) 

In addition, qualitative output corresponding judgment of the human eye was validated by employment 
of the questionnaire method among ten voluntary respondents who fulfilled the demographical criteria of: 
a) Age: between 20 and 24 years old. 

b) Health condition: excellent eyesight condition without the hindrance of colour blindness. 

This observing procedure is critical for triggering the correct assessment of brightness due to human 
eyes are sensitive to the light and colour space. Herewith, scaled responses were requested from each 
respondent upon appraising the brightness of virtually displayed images through sole judgment of their naked 
eyes. Engaged process concerning independent answering of the questionnaire has been further depicted in 
Figure 4. Each responder had one chance to decide for each of the 38 images. Each brightness assessment 
category; low, medium, and high was randomly arranged in the image sets. The setting in the room was 
controlled by using the same location for all of the respondents and the fluorescent lamp served as the light 
source. The display panel was calibrated by setting the brightness to zero. It is useful to avoid any concerns 
about the eye’s light adaptation that will decrease the accuracy of judgement. 


Identification Image Contrast 


r 


Luminance Arithmetic mean 


No 
Human visual 
observation 


Figure 3. Flowchart for image contrast measurement 
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Figure 4. Questionnaire responses from respondents 


2.3. Enhancement 

Figure 5 shows the entire flowchart used in this stage. Many researchers enhance the contrast 
images by applying the HE and CLAHE method without identifying the suitable images needed to apply the 
enhancement process. Therefore, this study proposed the pre-enhancement process that includes 2 steps: 
classification of the contrast scale category and contrast modification method. 


Pre-Enhancement 


a: Se 
Middle contrast 


" p i 


127-152 179-203 204-228 


Low contrast High contrast 


102-126 229-254 


Classification contrast scale category 


v 


Contrast modification method 
- Increase contrast 
- Decrease contrast 
(1096 - 5096) 


Contrast modification 


Figure 5. Block diagram for the enhancement process 


2.3.1. Identification of the contrast scale category 

The iris images were divided into three contrast levels: low, medium, and high. As demonstrated in 
Table 1, this classification was based on the Ansel Adam zone system [2]. Using the (1), the luminance 
value, Y; , was used to determine the contrast scale value and identify the grey name and contrast category. 


Table 1. Contrast category, scale value and gray name by Ansel Adam zone system 


Contrast category — Contrast scale value Gray name 
High 255 Pure white 
229—254 Bright white 
204-228 Gray 
Medium 179-203 Light gray 
153-178 Middle light gray 
127-152 Middle gray 
Low 102-126 Medium dark gray 
76-101 Very dark gray 
54-75 Dark black 
25-53 Near black 
0-24 Pure black 
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2.3.2. Contrast modification method 

The best contrast value for the human eye’s is around 60% to 70% from the maximum contrast scale 
value (255) [3], which is categorized as medium contrast. Therefore, this step aims to manipulate the contrast 
value up to the range of middle contrast scale value from 127-203. In the low contrast group, there are five 
grey names: pure dark, near black, dark black, very dark gray, and medium dark gray. All images in the low 
contrast category were divided into two groups based on scales: a) 0-101 and b) 102—126 for analysis. These 
two groups were formed since the results of the 0—101 contrast modification technique did not show any 
increments in contrast value. But, the results of images in the range 102—126 showed an improvement in 
contrast value. In the filter step, HE and CLAHE approaches were applied to increase the contrast value for 
images on the scale of 0—101. After the enhancement process, all the images on a scale of 0—101 need to 
undergo the contrast modification method, which involves increasing the contrast value 10%, 20%, 30%, 40%, 
and 50% from the original contrast. The red channel was utilised in this analysis in order to extract the data 
region of interest (ROD, which represents the colour of the blood vessels. After that, the contrast modification 
method was applied to the images in the middle and high contrast categories. The process is done by increasing 
the contrast by 10%, 20%, 30%, 40%, and 50% or decreasing the contrast by 10%, 20%, 30%, 40%, and 50%. 


2.4. Performance measurement 

The final step is the performance measurement between the enhancement and original images using 
mean square error (MSE) and peak signal to noise ratio (PSNR) shown in the (4) and (5), respectively. MSE 
is a technique to measure the mean value of error that occurs between the original image and the output image. 
While, PSNR represents the ratio between maximum value of contrast (R) and MSE of the output image. 


_ EMN[I Q4-N)-12 (M.N)]? 
= M*N 


MSE 


(4) 


Where M and N are denoted for number of rows and columns respectively. I, is the original image and I, is 
the enhanced image. 


2 
PSNR = 10 logro — (5) 


2.5. Validation 

The suggested algorithm had been verified using a contrast enhancement assessment dataset that 
was available to the public (CEED2016). The dataset was created to evaluate the efficacy of image 
enhancement [30]-[32]. The image in the dataset was divided into three categories according to its luminance 
value. Synthetic data was generated to increase the number of high contrast datasets by manipulating the 
contrast value of the middle contrast images. There are a total of 24 images for validation. 


3. RESULTS AND DISCUSSION 
3.1. Image contrast measurement results 

Table 2 shows the results of a similarity contrast measurement between brightness and arithmetic 
mean, which corresponded to the contrast scale value and gray name. The first image (43—264x357), for 
example, has a different contrast value between brightness and arithmetic mean, but the range in the middle 
gray is identical. Hence, the similarity status is similar. Figure 6 depicted a summary of similarity image 
contrast measurement between luminance and arithmetic mean for three groups. Group medium resolution 
yielded comparable results between similar and different. While for low and high resolution the status similar 
is lower than the status different. Finally, the similarity findings demonstrated that quantitative measurement 
is insufficient for presenting the characteristics of the gray image. 

Additionally, qualitative measurements as obtained from the visual observation of ten individual 
respondents on both luminance and arithmetic calculations have been outlined in Table 3. The qualitative 
results indicated that the human eye distinguishes colour could be easily detected using the luminance approach 
compared to the arithmetic mean. Therefore, the luminance technique was used for further processing. 


32. Enhancement results 

Table 4 shows the results of 38 iris images based on category, contrast scale, a gray name based on Ansel 
Adam zone system, resolution, luminance value before and after the pre-enhancement process. Red text indicates 
that it achieved the medium contrast category. Based on the luminance value after contrast modification method 
results in Table 4, images in the contrast scale 0—101 (pure black, near black, dark black, and very dark gray) did 
not achieve the optimum range scale. In addition, only three iris images that are very dark gray with contrast 
increment by 50% meet the medium category. Meanwhile, the resulting contrast modification method in medium 
dark gray (102—126) satisfies the medium contrast category after a contrast increment of 50%. 
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Figure 6. Summary of luminance and arithmetic mean contrast similarity measurement 


Table 2. Example result of contras measurement for low resolution image 


Low resolution Luminance (contrast value) Gray name Arithmetic (contrast value) Scale Status 
43-264x357 141 Middle gray 128 Middle gray Yes 
67—239x236 159 Middle light gray 128 Middle gray No 
15—208x208 89 Very dark gray 128 Middle gray No 
24-196x301 179 Light gray 128 Middle gray No 
82—261x400 117 Middle dark gray 128 Middle gray No 
94—183x275 73 Dark black 111 Middle dark gray No 

3—233x263 2 Middle gray 3 Middle gray Yes 
93—-232x360 78 Dark black 117 Middle dark gray Yes 

6—-325x350 58 Dark black 87 Very dark black No 
13-189x284 228 Medium dark gray 128 Medium dark gray Yes 
h4-185x273 60 Dark black 89 Very dark black No 
29-270x388 3 Pure black 6 Pure black Yes 
41-342x387 221 Gray 128 Middle gray No 
m4—220x229 61 Dark black 92 Very dark black No 
477—300x400 0 Pure black 0 Pure black Yes 
66—-342x400 1 Pure black 2 Pure black Yes 


Table 3. Respondents’ observation on contrast measurement 


Responder Low resolution image Medium resolution image High resolution image 
Luminance (%) Arithmetic (%) Luminance (%) Arithmetic (%) Luminance (%) Arithmetic (%) 
Responder 1 62.50 18.75 41.67 25.00 21.27 9.09 
Responder 2 56.25 18.75 41.67 25.00 18.18 0 
Responder 3 56.25 12.50 25.00 16.67 18.18 9.09 
Responder 4 56.25 18.75 33.33 16.67 9.09 9.09 
Responder 5 50.00 6.25 33.33 16.67 18.18 0 
Responder 6 56.25 18.75 41.67 25.00 18.18 0 
Responder 7 56.25 18.75 41.67 25.00 18.18 0 
Responder 8 56.25 18.75 41.67 25.00 9.09 9.09 
Responder 9 50.00 6.25 25.00 16.67 9.09 0 
Responder 10 43.75 6.25 33.33 16.67 27.27 9.09 
Average 54.38% 14.38% 35.83% 20.84% 17.27% 4.55% 


Previous results showed that the luminance value after the contrast modification method in the pure 
black, near black, dark black, and very dark gray (0—101) categories did not achieve the luminance value in 
the middle contrast category. Therefore, this study proposed two well known contrast enhancement methods 
which are HE and CLAHE. These methods applied only to all images in pure black, near black, dark black, 
and very dark gray categories. Table 5 shows the outcomes after using HE and CLAHE. 

Figure 7 shows a boxplot illustrating the further analysis of the results in Table 6. The boxplot 
displays the median score for CLAHE compares to the original images are slightly different. Therefore, it’s 
not shown a very significant image enhancement to the original images. Besides, the interquartile ranges 
between CLAHE and the original have a constant size. These indicate that the dispersion of luminance is 
low. Images with a luminance value greater than 100 have a high tendency to improve the contrast from low 
to the middle range. Meanwhile, the median score for HE has a significant output compare to the original 
images and the dispersion of luminance exceeding 100. Thus, these results showed that the HE method is a 
suitable technique compared to the CLAHE method for improving the luminance value in low contrast 
categories, particularly for scale values between 0-101. 
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Table 4. Result on enhancement process based on category 
Luminance value 
Category Scale Gray name Resolution Original Increase (%) Decrease (%) 

10 20 30 40 50 10 20 30 40 50 

High 255 Pure white - - - - - - - - - - - - 
229-254 Bright white 850x574 231 247 257 267 277 287 227 217 207 197 187 
204-228 Gray 819x1024 251 261 271 281 291 301 241 231 221 211 201 
189x284 228 238 248 258 268 278 218 208 198 188 178 
342x387 221 232 242 252 262 272 212 202 192 182 172 
Middle 179-203 Light gray 196x301 79 189 199 209 219 229 169 159 149 139 129 
153-178 Middle light gray 239x236 159 169 179 189 199 209 149 139 129 119 109 

127-152 Middle gray 264x357 41 152 162 172 182 192 132 122 112 102 92 

448x673 40 150 160 170 180 190 130 120 110 100 90 

600x772 137 146 156 166 176 186 126 116 106 96 86 

319x500 30 139 149 159 169 179 119 109 99 89 79 

Low 102-126 Medium dark gray 261x400 17 128 138 148 158 168 108 98 88 78 68 
640x640 09 119 129 139 149 159 99 89 79 69 59 

545x736 05 115 125 135 145 155 95 85 75 65 55 

080x 1080 03 113 123 133 143 153 93 83 73 63 53 

681x1000 03 113 123 133 143 153 93 83 73 63 53 

76-101 Very dark gray 492x749 96 106 116 126 136 146 86 76 66 56 46 

976x2704 81 91 101 111 121 131 71 61 51 41 31 

571x575 75 85 95 105 115 125 65 22 45 35 25 

54-75 Dark black 183x275 73 84 94 104 114 124 64 54 44 34 24 

220x229 61 71 81 91 101 111 51 41 31 21 11 

340x2016 61 71 81 91 101 111 51 41 31 21 11 

185x273 60 69 79 89 99 109 49 39 29 19 9 

325x350 58 68 78 88 98 108 48 38 28 18 8 

492x749 96 106 116 126 136 146 86 76 66 56 46 

976x2704 81 91 101 111 121 131 71 61 51 41 31 

25-53 Near black 393x521 33 42 52 62 72 82 22 12 2 -8 -18 
0-24 Pure black 809x1200 20 30 40 50 60 70 10 0 -10 -20 -30 

358x2048 18 28 38 48 58 68 8 -2 -12 -22 -32 

696x1024 10 21 31 41 51 61 1 -9 -19 -29 -39 

640x640 9 19 29 39 49 59 -1 -11 -21 -31 -41 
080x1080 ü 16 26 36 46 56 -4 -14 -24 -34 -44 

270x388 3 14 24 34 44 54 -6 -16 -26 -36 -46 

233x263 2 12 22 32 42 52 -8 -18 -28 -38 -48 

768x1024 2 12 22 32 42 52 -8 -18 -28 -38 -48 

342x400 1 11 21 31 41 51 -9 -19 -29 -39 -46 
300x400 0 10 20 30 40 50 -10 -20 -30 -40 -50 
450x600 0 10 20 30 40 50 -10 -20 -30 -40 -50 
1096x1280 0 10 20 30 40 50 -10 -20 -30 -40 -50 
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Table 5. Results after HE and CLAHE method 


Luminance value 


Value scale Resolution Original HE CLAHE 
76-101 492x749 96 98 101 
1976x2704 81 105 73 
54-75 571x575 75 115 148 
183x275 73 100 118 
220x229 61 227 67 
1340x2016 61 118 55 
185x273 60 117 54 
325x350 58 89 95 
25-53 393x521 33 152 37 
0-24 809x1200 20 0 13 
1358x2048 18 55 28 
696x1024 10 7 21 
640x640 9 0 6 
1080x1080 7 23 14 
270x388 3 2 14 
233x263 2 26 17 
768x1024 2 37 4 
342x400 1 24 5 
300x400 0 2 5 
450x600 0 12 3 
1096x1280 0 50 3 
Average 74.14 99.38 64.71 
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After that, the iris images after the HE method will go through the contrast modification method to 
improve the luminance value. Table 6 shows that the luminance value after contrast modification method in 
the (0-101) scale with HE method. From the results, images with the HE method after contrast increment by 
3096—5096 achieved the optimum contrast scale. HE method enhances the image contrast based on the entire 
image histogram with low computational. However, images in the 0—24 scale did not show any improvement 
for luminance value. The decrement contrast after HE method did not perform well at this stage due to the 
iris images in the (0—101) categories are not suitable for decrement contrast. 


Compare Luminance 


200 | 


100 | 


HE CLAHE ORIGINAL 


Figure 7. Luminance comparision between HE, CLAHE and original images 


Table 6. Results of enhancement process following application of HE technique 
Luminance value 

Scale Resolution . Increment (96) 
10 20 30 40 50 
76-101 492x749 96 98 108 118 128 138 148 
976x2704 81 105 115 125 135 145 155 
54-75 571x575 75 115 125 135 145 155 165 
183x275 73 100 110 120 130 140 150 
220x229 61 227 237 247 257 267 217 
1340x2016 61 118 128 138 148 158 168 
185x273 60 117 127 137 147 157 167 
325x350 58 89 99 109 119 129 139 
25-53 393x521 33 152 162 172 182 192 202 
0-24 809x1200 20 0 10 20 30 40 50 
1358x2048 18 55 65 75 85 95 105 
696x1024 10 T 17 27 37 47 57 


640x640 9 0 10 20 30 40 50 
1080x1080 7 23 33 43 53 63 73 
270x388 3 27 37 47 57 67 77 
233x263 2 26 36 46 56 66 76 
768x1024 2 37 47 57 67 TI 87 
342x400 1 24 34 44 54 64 74 
300x400 0 2 12 22 32 42 52 
450x600 0 12 22 32 42 52 62 
1096x1280 0 50 60 70 80 90 100 


3.3. Simulation results 

Finally, the performance measure is calculated based on the MSE and PSNR values. Figure 8 shows 
that the MSE and PSNR are low based on the recommendation method. The proposed method yields an average 
MSE 18.25 and an average PSNR of 28.87. Low MSE indicates that the image quality is good because an error 
in the image became decreasing while the higher PSNR value indicates the better quality of an image. Figure 9 
shows a visual output of the image enhancement. The original image is depicted in Figure 9(a), Figure 9(b) 
illustrates a contrast reduction of 1096, Figure 9(c) a contrast reduction of 2096, Figure 9(d) a contrast 
reduction of 30%, Figure 9(e) a contrast reduction of 40%, and Figure 9(f) a contrast reduction of 50%. 
On the other hand, Figure 9(g) presents a contrast increment of 10%, Figure 9(h) demonstrates a contrast 
increment of 20%, Figure 9(1) has a contrast increment of 30%, Figure 9(j) has a contrast increment of 40%, 
and Figure 9(k) depicts a contrast increment of 50%. 
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Figure 8. Average MSE and PSNR value 


(2) (h) Q) 


Q) (k) 


Figure 9. Iris image for pre-enhancement process with: (a) original image, (b) contrast decrement by 10%, 

(c) contrast decrement by 20%, (d) contrast decrement by 30%, (e) contrast decrement by 40%, (f) contrast 

decrement by 50% respectively, (g) contrast increment by 10%, (h) contrast increment by 20%, (i) contrast 
increment by 30%, (j) contrast increment by 40%, and (k) contrast increment by 50% respectively 
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3.4. Validation results 

Figure 10 illustrates the PSNR validation results for the enhancement of low and high contrast 
measurements. The graph has shown a consistence performance with a maximum PSNR is 74.68 and 
a minimum PSNR is 63.1. A greater PSNR indicates a good quality performance for the enhancement image. 
Meanwhile, the results obtained from the MSE are shown in Figure 11. The average value for low and high 
contrast being 0.012 and 0.15 respectively. It indicates that the proposed method produces high-quality images. 
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Figure 11. MSE validation 


4. CONCLUSION 

This study shows that many iris images have poor luminance value due to the low lighting during 
the capturing of iris images. The contrast adjustment in the pre-enhancement step is critical to avoid 
over-enhancement, which results in the loss of iris image information. In conclusion, a 50% increment of 
contrast is suitable for images with luminance values from the scale 102—126. The combination of HE 
method with a 50% contrast increment is appropriate for images scale 25—101. Meanwhile, for images in the 
middle contrast category, the original contrast is preserved. Finally, images in the high contrast category must 
decreasing contrast by 50%. The noise was minimised after the pre-enhancement procedure by contrast 
increment or decrement, resulting in low MSE and PSNR values. Besides that, the result shows that the 
proposed method was valid for all multiresolution images. In future, human assessment validation of the 
contrast category with scale and the gray name is needed to clarify the results. 
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